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Abstract: Let P = {hi, ...,hs}c ^Yi, . . . ,Yk], D > deg{h^) for l<i<s,a 
bounding the bit length of the coefficients of the hi's, and $ be a quantifier-free 
7^-formula defining a convex semi- algebraic set. We design an algorithm return- 
ing a rational point in S if and only if 5 n Q ^ 0. It requires o-C)(i)£)0(fc ) ^-^ 
operations. If a rational point is outputted its coordinates have bit length dom- 
inated by uD'^'^^ \ Using this result, we obtain a procedure deciding if a poly- 
nomial / e . . . , Xn] is a sum of squares of polynomials in Q[Xi, . . . , X„]. 
Denote by d the degree of /, t the maximum bit length of the coefficients in /, 
D = ("+'^) and k < D{D + 1) - This procedure requires r^i^) D^^''"^ 
bit operations and the coefficients of the outputted polynomials have bit length 
dominated by tD'^^'' K 
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Calcul de points rationnels dans des 
semi-algebriques convexes et decomposition en 
sommes de carres 

Resume : Soit V = {hi,...,hs} C Z[ri, . . . , Yfc], D > deg(/ii) pour 1 < 
i < s, a une borne sur la longueur binaire des coefScients des hi, et $ une V- 
formulc sans quantificatcurs dcfinissant un onsemble semi-algcbriquc convcxc. 
Nous dccrivons un algorithme qui retourne un point a coordonnees rationnelles 
dans S si et seulement si <S fl Q 0. Get algorithme est de complexite binaire 
(jO(i)]jO(k- ) gj point rationnel est renvoye, ses coordonnees sont de longueur 
binaires dominees par crD^C^ ). On deduit de ce resultat une procedure qui 
decide si un polynfime / G Z[Xi, . . . , X„] est une somme de carres de polynSmes 
dans Q[Xi. .... X,,]. Soit d \c, dcgre de /, r le maximum des longueurs binaires 
des coefficients de /, D = et fc < D{D + 1) - {"+^') . Cette procedure est 

de complexite binaire t^C^^D^^'' ^ et les coefficients des polynfimes obtenus en 
sortie ont une longueur binaire dominee par tD'^'^^ \ 

Mots-cles : sommes dc carres a coefficients rationnels, programmation semi- 
definie positive, ensembles semi-algebraiques convexes, complexite. 
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1 Introduction 

Motivation and problem statement. Suppose / £ M[a;i, . . . , a;„], then / is 
a sum of squares (SOS) in R[xi, . . . , x„] if and only if it can be written in the 
form 

/^v^-M-v, (1) 

in whic h v is a column vector of mono mials and M is a real positive sera i defini te 
matrix ( Powers and Wormannl . 19981 . Theorem 1) (see also Choi et al. (1995)). 



M is also called a Gram matrix for /. If M has rational entries, then / is a sum 
of squares in , . . . , x„] . 

PROBLEM 1.1 (Sturmfels). If f ^ Q[a;i, . . . , a;„] is a sum of squares in 
R[a;i, . . . , Xn], then is f also a sum of squares in . . . , Xn] ? 

It has been pointed out that if there is an invertible Gram matrix for /, then 
there is a Gram matrix for / with rational entries fjHillari . .2009. Theorem 1.2). 
Furthermore, if / G (Q)[a;i, . . . , Xn] is a sum of m squares m K[xi, . . . , x„], where 
K \s, & totally real number field with Galois cl osure L, the n / is also a sum of 
4m • 2[^^01+i([^^^l+^) squares in Q[xi, . . . , a;„] (Hillarl. l2009l . Theorem 1.4). It is 
i ntere sting to see that the number of squares can be reduced to m (see KaltoferJ 
()2009l )V 



Although no example is known of a rational polynomial having only irra- 
tional sum of squares, a complete answer to Question [TTT] is not known. This is 
the main motivation for us to design an algorithm to check whether a rational 
polynomial having a rational sum of squares decomposition and give the rational 
SOS representation if it does exist. By reducing this problem to semi-definite 
programming, this can be done by designing a n algorithm checking if a conv ex 



semi-algebraic set contains rational points (see lPowers and Wormannl ||199 



Main result. We propose an algorithm which decides if a convex semi-algebraic 
set S <zM.^ contains rational points (i.e. points with coordinates in Q^). In the 
case where S n Q*^ is non-empty, a rational point in S is computed. 

The semi-algebraic set S is given as the solution set of a polynomial system 
of non-strict inequalities with integer coefficients. Arithmetic operations, sign 
evaluations and comparisons of two integers/rationals can be done in polynomial 
time of the maximum bit length of the considered integers/rationals. 

We bound the number of bit operations that the algorithm performs with 
respect to the number of polynomials, their degrees and the maximum bit length 
of their input coefficients; we also give upper bounds on the bit length of the co- 
ordinates of the outputted rational point if this situation occurs. More precisely, 
the main result is as follows. 

Theorem 1.1 Consider a set of polynomials V — {hi, . . . , hs} C Z[Yi, . . . , Yfe], 
and a quantifier-free V -formula '^{Yi, . . . , Yk) and let D he an integer such that 
deg{hi) < D for 1 < i < s and a the maximum bit length of the coefficients 
of the hi's. Let S C M'"' be the convex semi- algebraic set defined by <&. There 
exists an algorithm which decides ifSCiQ'' is non-empty within a'^'^^^sD)'^^'' ^ 
bit operations. In case of non-emptiness, it returns an element ofSnQ'^ whose 
coordinates have bit length dominated by aD'^'^^ ^ . 
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We use a procedure due to Basu et al. 1 1996l l performing quantifier elimina- 
tion over the reals in order to deduce from Theorem [TTT] the following result. 



Corollary 1.2 Let S C M.'' be a convex set defined by 
S = {Y eR'' : (QiXW e M"i) • • • e I 



P{Y,X 



[1] 



vrith quantifiers Qi G {3,V}, where is a set of ni variables, P is a Boolean 
function of s atomic predicates 

5(y,xw,...,xM)A,o 

where S {>, <, =} (for i = 1, . . . ,s) and the gi 's are polynomials of degree 
D with integer coefficients of binary size at most a. There exists an algorithm 
which decides ifSdQ'' is non-empty within cr*^^^^ (sZ?)'-*^*' ^i=i'^i) jjH operations. 
In case of non-emptiness, it returns an element ofSOQ'' whose coordinates have 
bit length dominated by aD'^'^''^^'i'=^'^'\ 

The proof of the above results is based on quantitative and algorithmic re- 
sults for computing sampHng points in semi-algebraic sets and quantifier ehmi- 
nation over the reals. 

It is well-known that deciding if a given polynomial / G . . . , X„] of 

degree d whose coefficients have bit length dominated by t is a sum of squares of 
polynomials in Q[Xi, . . . , X„] can be reduced to a linear matrix inequali t y whi ch 



Powers and Wormann 1 199. 



defines a convex semi-algebraic set (see e.g. 
Applying Theorem II. !( we show that there exists an algorithm deciding if such 
an SOS decomposition exists over the rationals and that the coefficients of the 
polynomials in the decomposition have bit length dominated by tD'^^'' ^ with 
and k < D{D -|- 1) — ("^^''). Moreover, such a decomposition can 



D = 



ln+d\ 



be found within t'^W ]jO(k ) operations. 



Prior works. Khachiyan and Porkolab extended the well-known result of 
Lenstral 1 1983l l on the polynomial-time solvability of linear integer program- 
ming in fixed dimensi on to semidefinite integer programmi ng. The following 
proposition is given in 



m to semidehnite mteger programmi i 
Khachiyan an d PorkolabI l|l997l . l20Q0f ) 



Proposition 1.3 Let S C M. be a convex set defined as in Corollary 
There exists an algorithm for solving the problem minlY^jy = (Yi, . . . ,Yk) G 
Sf^Z''} in time i'^'^^^sD)^'^'' )njLiO(ni) ^^^g of non-empty, then the min- 
imization problem has an optimal solution whose bit length is dominated by 

Their algorithm was further improved by Heinz for the case of convex min- 
imizat i on w here the feasible region is described by quasiconvex polynomials 
Heinj l|2005l l. 

Although we can apply Proposition 11.31 directly to certify that a given poly- 
nomial with integer coefficients to be non-negative for all real values of the 
variables by computing a sum of squares in Z[a;i, . . . , the nonnegativity 
of a polynomial can be certified if it can be written as a sum of squares of 

polynomials in 0[a:i x^]. Some hybrid symbolic-numeric algorithrn s have 

been given in lPevrl and Parrild (|2007l . l2008l l: lKaltofen et all l|2008l . l2009[ ) which 
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turn a numerical sum of squares representation of a positive polynomial into 
an exact rational identity. However, it is well known that there are plenty of 
polynomials which are nonnegative but can not be written as sums of squares of 
polynomials, for example, the famous Motzkin polynomial. This also impel us 
to study Khachiyan and Porkolab's approach. It turns out that by focusing on 
rational numbers instead of integers, we can design an exact algorithm which 
decide whether a given polynomial can be written as an SOS over the rationals 
and give the rational SOS decomposition if it exists. 



Structure of the paper. Section [2] is devoted to recall the quantitative 
and algorithmic results on computing sampling points in semi-algebraic sets 
and quanti fi er eli mination over the reals. Most of these results are proved in 
Basu et al. ( 1996f ). Section [3] is devoted to prove the correctness of the algo- 



rithm on which Theorem 1 1.1 1 and Corollary 11.21 rely. The complexity analysis is 
done in Section HI In Section \5\ we apply Theorem 11.11 to prove the announced 
bounds on the bit length of the rational coefficients of the decomposition into 
sums of squares of a given polynomial with integer coefficients. 



Acknwledgments. This work is supported by the EXACTA grant of National 
Science Foundation of China (NSFC) and the French National Research Agency 
(ANR). The authors thank INRIA, KLMM and the Academy of Mathematics 
and System Sciences for their support. 



2 Preliminaries 

The algorithm on which Theorem 11.11 relies and its complexity analysis are 
based on algorithmic and quantitative results on computing sampling points in 
semi-algebraic sets and quantifier elimination over the reals. 



2.1 Computing points in semi-algebraic sets 

Consider a set of polynomials V = {hi, . . . , hj} C Z[Yi, . . . , Y^], and a quantifier- 
free "P-formula $(^1, . . . , Yk) (i.e. a quantifier-free formula whose atoms is one 
of /i = 0, ft- 7^ 0, /i > 0, < for /i G 7^). Let D be an integer such that 
deg(/ii) < D for 1 < « < J and £ the maximum bit length of the coefficients of 
the hi's. We denote by iS C K'^ the semi-algebraic set defined by ■ • • , i^A 



k ■ 



A function RealizableSignConditions computing a set of algebraic points hav- 
ing a non-empty intersection with each connected componen t of semi-algebraic 
sets defined by sig n conditions satisfi ed by V is given in I Basu et al . 199(1 



Section 3) (see also ( Basu et al. . 2Q06I . Chapter 5)). From this, a function Sam 



plingPoints computing a set of algebraic points having a non-empty intersection 
with each connected component of S is obtained. These algebraic points are 
encoded by 



a rational parametrization 
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where G, Go, . . . ,Gk are polynomials in Z[T] such that deg(gcci(G', Go)) = 
and 

for 1 < i < fc, -1 < deg(G,) < deg(G) - 1 and < deg(Go) < deg(G) - 1; 

the rational parametrization is given by the list Q — (G, Go, Gi, . . . , G^); 
the degree of G is called degree of the rational parametrization and Z{Q) (Z 
C'^ denotes the set of complex points encoded by G', 

and a list T of the Thom-encodings of the real roots i9 of G such that 
$ (£im £A^] is true 



The bit complexity of SamplingPoints is iJ'^'^^D'^^'^^ and the output is such that 
deg(G) = 0{D)^ and the bit length of the coefficients of G, Gq, Gi, . . . , Gfc is 
dominated by W°'^''\ 

Factorizing over Q a univariate polynomial ft. € Q[r] of degree 5 with rational 
coeff icients of max i mum b it length £ can be done in £0(i)<50(i) bit -operations 
(see iLenstra et all l|l982[ ): Ivan Hoeii and NovocinI l|2007[ ): ISchonhagel ||l984)). 



Given a root ■() of h, the mi nimal polynomia l of d has coefficients of bit length 
dominated by t + 0{5) fsee lMignott^ l|l982l ll. 



Consider now a root -d oi G and its minimal polynomial g. Since G and Go 
are co-prime, one can compute G(^^ mod g to obtain a rational parametrization 
(5, go, • ■ ■ , 5fc) with integer coefficients of bit length dominated by iD'^^''^ and 

for 1 < i < fc, — 1 < deg{gi) < deg{g) — 1 and < deg(go) < deg(,g) — 1 

within a bit-complexity |?C)(i)£)0(fc) ^ This implies the following result. 

Proposition 2.1 There exists a function SemiAlgebraicSolve which takes as 
input the system <i>(li, . . . , Y/c) and computes a rational parametrization Q = 
(G, Go, Gi, . . . , Gfc) and a list T of Thom-encodings such that G is irreducible 

overQ, andT contains the encodings of the real roots d ofG such that ■ • ■ , §^py 

S. The hit length of the coefficients of G, Go, Gi, . . . , Gfc is dominated by iD'^'^^^ 
and deg(G) = 0{D)^ . Moreover, SemiAlgebraicSolve requires jfe+i^o(fc) 
hit operations. 

Remark 2.2 Since G and Gq are co-prime, one can compute Gq^^ mod G in 
polynomial time, and the binary length of its rational coefficients can be bounded 
via subresultants, we can assume, without loss of generality, that the rational 
parametrization has a constant denominator: 

Y - i(Gi(^9), G2 . . . , Gkm e S, G(^) - 0, (2) 

where the bit length of q and the coefficients 0/ G, Gi, . . . , Gfc are dominated by 

The above discussion leads also to the following result. 

Proposition 2.3 Let Q,T be the output of SemiAlgebraicSolve(<i>), S be the 

degree of G, and i be the maximum bit length of the coefficients of the poly- 
nomials in Q U V. There exists a function RationalZeroDimSolve which takes 



G 
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as input Q and $ and returns a rational point y G Z{Q) if and only if y ^ 
S n Z(Q) n Q*^, else it returns an empty list. The coordinates of these ratio- 
nal points have bit length dominated by £5'^^-^'> and computations are performed 

within 0(fc)0(J)^o(i)<5O(i)("+^)°^'^ bit operations. 

Remark 2.4 According to Proposition [2J\ i/ie /i/nciion SemiAlgebraicSolve com- 
putes a rational parametrization Q ~ (G, Go, Gi, . . . , Gfe) such that G is irre- 
ducible over Q. Therefore a rational point y G Z{Q) if and only j/deg(G) = 1. 
In order to check whether y G S, we only need to evaluate the formula $ at y. 



The following result is a restatement of i Basu et al . 199(1 Theorem 4.1.2 



and allows us to bound the bit length of rational points in non-empty semi- 
algebraic sets defined by strict polynomial inequalities. 

Proposition 2.5 Let S' C be a semi- algebraic set defined by a quantifier- 
free P -formula whose atoms are strict inequalities. Then S' contains a rational 
point whose coordinates have bit length dominated by iD'-'^'^^ . 



The proof of the above result (see i Basu et al.l . Il99d Proof of Theorem 



4.1.2 pp. 1032)) is based on the routine RealizableSignConditions and the iso- 
l ation of real root s of univariate polynomials with rational coefficients (see e.g. 
( Basu et aD . l2006l . Chapter 10)). We denote by RationalOpenSemiAlgebraicSolve 



a function taking as input the "P-formula $ and which returns a rational point 
in S if and only if there exists a non-empty semi-algebraic set S' defined by a 
quantifier-free T'-formula whose ato ms are strict in equalities such that S' C S. 
The result below is not stated in iBasu et al. ( 1996h but is an immediate conse- 
quence of this proof. 

Corollary 2.6 Suppose that there exists a quantifier-free V -formula whose atoms 
are strict inequalities defining a non-empty semi-algebraic set S' C S. There 
exists an algorithm computing a rational point in S if and only if S ^ 9. It 
requires ^^(i) j*;+i£)0(fc) j^-^ operations and if a rational point is outputted, its 
coordinates have bit length dominated by iD'^^'^K 

2.2 Quantifier elimination over the reals 

We consider now a first-order formula F over the reals 

(Qi^I^l gM"^)---(Qc.x["1 £R"-) P{Y,X^^\...,X^'^^) 

where 

• Y = (Yi, . . . , Yfc) is the vector of free variables; 

• each Qi {i = 1, ... ,Lu) is one of the quantifiers 3 or V; 

is a Boolean function of s atomic predicates 

5(y,xw,...,xM)A^o 

where G {>, <, —} (for i = 1, . . . , s) and the gi^s are polynomials of 
degree D with integer coefficients of binary size at most £. 
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The following result on quantifier elimination is a restatement of (jBasu et al 

19961 . Theorem 1.3.1). 



Theorem 2.7 There exists a quantified-free formula 5" 

i=ij=i 

(where hij G "^[Yi, . . . , Yk] and e >} ) which is equivalent to F and such 
that 

m I < s(fe+i)nr=i(»i+i)i:)(fc+i)nr=iO(n.)^ 
. J, < snr=i(n,+i)£,n-^,oK)^ 

• deg{h,j) < i:)nr=iO(".)^ 

• the bit length of the coefficients of the polynomials hij is dominated by 
£23(fc+i)nr=iO(ni)_ 

The above transformation requires £s('=+i)nfai(ni+i)£)(fc+i)n^^iO(ni) j^-^ gpg^g^. 
tions. 

In the sequel, we denote by QuantifierElimination a function that takes F as 
input and returns a Hst [^'i, . . . , i] where the ^[s are the conjunctions 

J. 

/\{h,, A,,0). 



3 Algorithm and correctness 

3.1 Description of the algorithm 

We use the following functions: 

• Substitute which takes as input a variable Y^ G {Yi, . . . , Yfe}, a polynomial 
h G Q[yi, . . . , Yk] and a Boolean formula F and which returns a formula 
F obtained by substituting Yr by h in F. 

• RemoveDenominators which takes as input a formula F and returns a for- 
mula F obtained by multiplying the polynomials in F by the absolute 
value of the 1cm of the denominators of their coefficients. 

Consider now a rational parametrization Q = (G, Gq, Gi, . . . , Gfc, Gk+i) C 
Z[T]'=+3 with 5 = deg(G). For < i < (5 - 1, denote by a, G ll" the vector 
of integers whose j-th coordinate is the coefiicient of in Gj . Similarly, for 
Q < i < 5 — l,hi denotes the coefiicient of T' in Gk+i- We use in the sequel a 
function GenerateVectors that takes as input a rational parametrization Q. This 
function returns the set Hst of couples (a^jb^) for < i < 5 — 1. 

As in the previous section, consider now a set of polynomials V — {hi, . . . , hs\ C 
Z[Yi, . . . , Yfe], and a quantifier-free "P-formula $(Yi, . . . , Y^) and let D be an in- 
teger such that deg{hi) < D for 1 < i < s and a the maximum bit length of the 
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coefficients of the hi's. We denote by <S C ffi*^ the semi-algebraic set defined by 
$(Yi, . . . , Yfe) which is supposed to be convex. 

The routine FindRationalPoints below takes as input the formula ^{Yi, . . . , Yk) 
defining <S C M*^ and the list of variables [Vi , . . . , Yfe] . 

FindRationalPoints($, [Yi, Yk\). 

1. Let L = RationalOpenSemiAlgebraicSolve(Open((l>)) 

2. If X is not empty then return L 

3. Let g,T = SemiAlgebraicSolve($) 

4. If T is empty then return [] 

5. Let L = RationalZeroDimSolve(S, $) 

6. If Z/ is not empty or A; = 1 then return L 

7. Else 

(a) Let Ai, . . . , Ak, B be free variables and G be the formula 

Vy e M*^ + • • • + A| > A V AiYi + ■■■ + AkYk = B) 

(b) Let [^'i, . . . , */] = QuantifierElimination(0) and i = 1 

(c) While i < / do 

i. g,T = SemiAlgebraicSolve(4'i) and (G, Go, Gi, . . . , G^, Gfe+i) = 

g 

ii. If T is empty i = i + \ else break. 

(d) Let C = GenerateVectors(G, Go, Gi, . . . , G^, G^+i) 

(e) Let a = (ai, . . . , au) 7^ (0, . . . , 0) and S Z such that (a, b) € C 

(f) Let r — max(i, 1 < i < k and 7^ 0) 

(g) Let h = b- 

(h) Let $' = RemoveDenominators(Substitute(yr, ^, 'f)) 

(i) Let L= FindRationalPoints($',[yL,...,yr-i,i^r+i,---,ife]) 
(j) If L is not empty, 

i. Let (g'l, . . . , 5^-1, 9r+i, ■ ■ ■ ,<lk) be its element; 

ii. Let Qr = Evaluate({yi = Qi, 'i- < i < k, j ^ r}, h) 

iii. if $((7i, . . . ,qr-i,qr,qT+i, ■■■,qk) is true, return [{qi, . . . , g^-i, 9r, ^r+i, • • • ,9fe)] 
else return []. 

(k) Else return []. 

Proposition 3.1 The algorithm FindRationalPoints returns a list containing a 
rational point if and only if S fl Q*^ is non-empty, else it returns an empty list. 

The next paragraph is devoted to prove this proposition. 
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Remark 3.2 LetS cR'' be a convex set defined by 

S^{Y eR'' :M'^(giXW G M"0 • • • (Qo;^''^' £ K"-) p(r, X™, . . . , xi'^i)} 

vnth quantifiers Qi G {3,V}, where is a set of Ui variables, P is a Boolean 
function of s atomic predicates 

.9(y,xW,...,xM)A,0 

where G {>,<,=} (for i — 1, . . . ,s). 

Denote by the quantified formula defining S and by [^'i, . . . , ^ i] the output 
o/ QuantifierElimination(6). Running FindRationalPoints on the ^i's allows to 
decide the existence of rational points in S. This proves a part of Corollary \1.2[ 

3.2 Proof of correctness 



In th e sequel, we denote by closzar('5) its Zariski-closure. Following ijBochnak et al 



19981 . Definition 2.8.1 and Proposition 2.8.2 pp. 50), we define the dimension of 



S as the Krull dimension of the ideal associated to closzar(>5). By convention, 
the dimension of the empty set is —1. 

We reuse the notations introduced in the description of FindRationalPoints. 
The proof is done by induction on k. Before investigating the case fc = 1, we 
recall some elementary facts. 



Preliminaries. 

We start with a lemma. 

Lemma 3.3 Let A C M*'' be a semi-algebraic set defined by a quantifier-free 
V -formula. If dm\{A) = k there exists y G M*^ such that for all h h{y) > 
or h{y) < 0. 

Proof. Suppose that for all y € A, there exists h E V such that h{y) = 0. 
Then, A is contained in the union Ti, of the hypersurfaces defined by /i = for 
h V. Consequently, dim(^) < dim(7i) < k, which contradicts dim(^) = fc. □ 
The following lemma recalls an elementary property of convex semi-algebraic 
sets of dimension 0. 

Lemma 3.4 Let A C M.'^ be a convex semi- algebraic set. If dim(A) = Q, then 
A is reduced to a single point. 

Proof. If there exist two distinct points yi, 2/2 in A, the set B = {tyi + (1 — 
t)y2,t G [0, 1]} is contained in A. This implies that closzar(^) C closzar(^) and 
consequently dim(i3) < dim(A). Since closzar(-B) is the line containing j/i and 
1/2, dim(i?) = 1 and dim(A) > 1 which contradicts the assumption dim(yl) = 0. 
Our claim follows. □ 



Correctness when k = 1. 

Lemma 3.5 Suppose that k = 1. Then Steps (Qlt^j return a rational point in 
S if and only if S O Q'^ ^ else an empty list is returned. 

Proof. If fc = 1, the dimension of S is either 1, —1 or 0. 
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1. Suppose that S has dimension 1. From Lemma f3.3[ there exists a non- 
empty semi-algebraic set S' C S defined by a quantifier-free P-formula 
whose atoms are strict inequaUties. Thus <S' contains a rational point. 
From Corollary 12.61 such a rational point in S is outputted at Step |(T]). 

2. Suppose that S has dimension —1 (i.e. S is empty). From Proposition 
12.11 the Hst of Thom-encodings outputted at Step ^ is empty and the 
empty list is returned at Step 

3. Suppose that S has dimension 0. From Lemma [3^ 5 is a single point 
contained in Z{Q). From Proposition 12. 3[ this point is outputted at Step 
(lU if and only if it is a rational point; else the empty list is outputted. 

□ 

The case k> 1. 

Our induction assumption is that, given a quantifier-free "P'-formula (with 
V C . . . , Yfe_i]) defining a convex semi-algebraic set S' C M*^"^, FindRa- 

tionalPoints returns a Hst containing a rational point if and only if <S' fl Q'^^^ is 
non-empty, else it returns an empty list. 

Lemma 3.6 Suppose that < dim(5) < k. There exists (ai, . . . , a^) G K'^ and 
6 G M such that (ai, . . . , Ofe) ^ (0, . . . , 0) and 

V(yi,...,?;fe) gR*-' {yi,...,yk)eS=>aiyi + --- + akyk^b. (3) 

Proof. It is sufficient to prove that closzar('5) is an affine subspace over M: 
in this case, there exists a real affine hyperplane H (defined by X]i=i ^i^i = ^ 
for (ai, . . . , flfe) G R*^ \ (0, . . . , 0) and & G R) such that S C closzar(<S) C H. 

We prove below that closzai ('5) flR'^ is an affine subspace which implies that 
closzar('5) is an affine subspace. 

From Lemma [374^ if dim(5) — then 5 is a single point; thus the conclusion 
follows immediately. 

We suppose now that dim (5) > 0; hence S is not empty and contains in- 
finitely many points. Consider yo G S. Given y G R'^ \ {yo}, we denote by 
■^yo-.v the real line containing y and j/o and by -ffyo,y C R*^ the real affine 

hyperplane which is orthogonal to Lyp.a and which contains j/o- 

Since S is convex, for all y G 5\{yo}5 SCiLy^y^ ^ 0. We consider the set Uy^ — 
r\y^s\{yo} ^yo,y'^ note that Uy^ is an affine subspace since it is the intersection 
of affine subspaces. We claim that the orthogonal of Uy„ is closzar('5) n R'^. 

We first prove that S is contained in the orthogonal of Uy^ which impHes that 
closzar('5)nR'^ is Contained in the orthogonal oiUy„. By definition oiUyg, for all 
u G Uyg and all y G iS\ {j/o}, the inner product of you and yo, y is zero. We prove 
now that the orthogonal ofUyg is contained in closzar('5)nR'^'. By definition, the 
orthogonal oiUyg is the set of lines Ly^y^^ for y ^ S \ {yo}. Thus, it is sufficient 
to prove that for all y G 5 \ {yo}, Ly^y^ is contained in closzar('5) Pi R*^. For all 
y G 5\{j/o}, S{^Ly^yg ^ because 5 is convex. Moreover, closzar('5nij^^j^o)nR'^ 
is Ly^yg. Since SOLy^yg C S, Ly^yg is contained in closzar('5)nR'^. Our assertion 
follows. □ 

Suppose that dim(5) = k. Then, by Lemma [3731 S flQ'^ is not empty and a 
rational point is outputted at Step iJH) by Corollary [2761 Suppose now that S is 
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empty. Then, an empty list is returned at Step We suppose now that S is 
not empty and that no rational point is outputted at Step 1^. Hence, we enter 
at Step O. 

Remark that the formula 8 (Step l(7a|) l defines the semi-algebraic set A C 
M*^ X M such that (oi, . . . , Ofc, 5) G ^ if and only if (ai, . . . , ak) 7^ (0, . . . , 0) and 

V(?/i, . . . , yfc) e M'' (2/1, . . . , yfc) e 5 =^ aiyi H h akVk = b. 

Thus, the quantifier-free formula Vf=i (Step ifZbl) ) defines ^. Note that by 
Lemma 13.61 is not empty. Hence, the loop at Step ((7cl) ends by finding 
a rational parametrization G = {G,Gq,Gi, . . . ,Gk,Gk+i) (computed at Step 
(7(c)i|) which encodes some points in A. 

From the specification of SemiAlgebraicSolve, G is irreducible over Q. Let 
a = (ai, . . . , Ofc) G K'' and 6 G M such that (a, 6) G Ar\Z{Q). Then, there exists 
a real root d oi G such that 

, . dcg(G)-l , , 

where the couples (a^jb^) (ill' xT, are those returned by GenerateVectors (Step 
(fTd]) ). Since gcd(G'o,G) = 1, Go[-d) ^ 0. Moreover, (a, 6) G implies a ^ 
(0, . . . , 0). Note also that (a, 6) G implies that for all A G M*, (Aa, A6) G A 
since for all (yi, . . . , yu) G 5 and A G K* 

aiVi H h auVk = b <f=^ A(aiyi H h 0^?/^) = Xb 

This proves that 

(a*, 6*) = {GoWa, Go^b) e A and (a^, . . . , 4) ^ (0, . . . , 0). 

Thus, there exists i such that 0, which impHes that Step l(7e|) never fails. 
To end the proof of correctness, we distinguish the case where S CiQ'' is empty 
or not. 



The non-empty case. We suppose first that SDQ'^ is non-empty; let {yi, . . . , yt) G 
iSnQ'^. Using (U, the linear relation a\yi + - ■ ■ + a^yk — b* implies the algebraic 
relation of degree deg(G) — 1: 

deg(G)-l fc 

^ ^*(^a,,,y,-bO = 0, (5) 

i=0 j=l 

where anj is the j-th coordinate of a^. Since G is irreducible, it is the minimal 
polynomial of hence i9 is an algebraic number of degree deg(G). Thus, ^ is 
equivalent to 

fc 

VO < z < dcg(G) - 1, Yl ^'iVj = 
i=i 

We previously proved that there exists i such that a^ 7^ 0. We let a = 
(oi, . . . , flfc) G Z'^' \ (0, . . . , 0) and 5 G Z be respectively the vector with inte- 
ger coordinates and the integer chosen in G (Step l(7e|) ). We have just proved 
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that S n Q'^ is contained in the intersection of S and of the afHne hyperplane 
H defined by aiFi + • ■ • + UkYk = b. Note also that S H H is convex since S is 
convex and H is an affine hyperplane. 

Consider the projection TTr : (yi, . . . , j/fc) £ M'' ^ (yi, . . . , j/^-i, y^+i, ■ ■ ■ ,yk) <^ 
for the integer r computed at Step ([TJ). It is clear that the formula 
computed at Step l|7hp defines the semi-algebraic set TiriS n H) C R'^"^. Since 
5ni/ is convex, iTr{SC\H) is convex. Thus, the call to FindRationalPoints (Step 
(fTH) ) with inputs and [Yi, . . . , y^-i, i^r+i, ■ • ■ , ^fe] is valid. From the induction 
assumption, it returns a rational point in 77^(5 fl H) if and only if ■Kr{S fl H) 
has a non-empty intersection with Q'^~^. 

Since 5 n Q*^ (which is supposed to be non-empty) is contained in 5 fl iJ, 
-Kris n H) contains rational points. Thus, the list L (Step |[7l])) contains a 
rational point qfe_i = {qi, . . . , g^-i, 9r+i, ■ ■ ■ ,qk) G 7rr(5 n i/). This implies 
that 7r~^(qfc_i) fl H has a non-empty intersection with S H H. Remark that 
7r~^(q/c_i) ni/ is the rational point q = (qi, . . . , g^-i, ^r, 9r+i, ■ ■ ■ ,qk) where qr 



is computed at Step (7(j)ii|. It belongs to S since 7r~^(qfe_i)ni? and SnH have 
a non-empty intersection. Thus, $(gi, . . . , (/r-i, 9r, ^r+i, ■ • ■ , is true and q is 
returned by FindRationalPoints. 

The empty case. Suppose now that S DQ^ is empty. As above H denotes 
the affine hyperplane defined by aiYi -I- • • • -I- OfcYfe = b where (ai, . . . , at) S Z*"' 
and 5 e Z are chosen at Step l(7el). Using the above argumentation, nr{S n H) 
is convex and the formula $' (Step ((7h|) ) defines nr{S fl if). Thus, the call 
to FindRationalPoints (Step ^) with inputs and [Yi, . . . , Fr-i, Y^+i, . . . ,yfc] 
is valid. Suppose that Trr{S n H) does not contain rational points. Then, by 
the induction assumption, L is empty and the empty list is returned (Step 
(T^) which is the expected output since we have supposed 5 n Q*^ = 0. Else, 
L contains a rational point (gi, . . . , qr-i, qr+i, ■ ■ ■ , qu)- Consider the rational 



point (gi, . . . , qr-i,qr, qr+i, ■ ■ ■ ,qk) (where qr is computed at Step (7(j)ii|). It 
can not belong to S since we have supposed S n Q'' is empty. Consequently, 
^{qi, . . . , Qr-ij qr, qr+i, ■ ■ ■ , qk) is false and the empty list is returned. 

4 Complexity 

We analyze now the bit complexity of FindRationalPoints. 

Proposition 4.1 Consider a set of polynomials V = {/ii, . . . , h^} C Z[Yi, . . . , Yk] 
and a quantifier-free V -formula ^{Yi, . . . , Yk) and let D be an integer such that 
deg(/ii) < D for 1 < i < s and a the maximum bit length of the coefficients of 
the hr's. Then, FindRationalPoints($, [Yi, . . . , Yfc]) requires a°^^\sD)°^'''") bit 
operations. Moreover, if it outputs a rational point, its coordinates have bit 
length dominated by aD'^^'' K 

Remark 4.2 Let S C M.''' be a convex set defined by 

S = {Y eR'' -.R'^iQiX^^^ eM"i)---(Q^x['^l eR"-)P(r,xW,...,x["l)} 

with quantifiers Qi £ {3,V}, where X'*! is a set of Hi variables, P is a Boolean 
function of s atomic predicates 

5(Y,XW,...,XH)A,0 
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where Ai G {>, <, —} (for i = 1, . . . ,s) and the gi 's are polynomials of degree 
D with integer coefficients of binary size at most a. Denote by 8 the quan- 
tified formula defining S. By Theorem \2.7[ QuantifierElimination(6) requires 
crs('=+i)nr=i(«.+i)i:»(fc+i)nr=iO(«,) operations. 

It outputs a list of conjunctions $i, . . . , $/ withi < s(fc+i)n"=i("i+i)£)(fc+i)n^=iO(n,)^ 
and for 1 < i < I , is a conjunction of Ji < s'^i=i^"'+^^_D'^>=i°'^"') atomic 
predicates HAO with h G Z[Yi, . . . , Yfc], A G >} and deg{h) < Dnr=iO(«.) 
and the bit length of the coefficients of the polynomials hij is dominated by 
(jX}('=+i)n"=iO(n,) Thus, the cost of running FindRationalPoints on all the ^i's 
requires a'^'^^^sD)'^'^^ n.^^ni) operations. In case of non-emptiness of Sf]^^ , 
it returns an element of S H Q'^ whose coordinates have bit length dominated by 
(j£)0{k n^^in,) rpj^^g g^^^ prove Corollaru \l.2l 

We start with a lemma. 

Lemma 4.3 Steps mUB) o/FindRationalPoints(<I>) perform within a'^^^^ s'^^'^ D'^^'''> 
bit operations. If a rational point is returned at Step (0j or Step l^jj, its coordi- 
nates have bit length dominated by aD'^'^^^ . 

Proof. The result is a direct consequence of the results stated at Section [H 

1. From Corollary 12. 6t Step |[T]) is performed within us'^+^D'^^'^' bit opera- 
tions and if a rational point is outputted at Step iJH) , its coordinates have 
bit length dominated by crD'-'^*'^ . 

2. From Proposition l2.1l Steps ^ and (0]) are performed within 0-0(1) gfe+i£)0(fe) 
bit operations. 

3. From Proposition [231 Step ([5]) requires a'^'^^^ D'^'^^^ bit operations. More- 
over, if a rational point is outputted at Step |[6]), its coordinates have bit 
length dominated by aD'^^''\ 

□ 

We prove now the following result. 

Lemma 4.4 1. Steps f 7aH7/tp require a'-' (i)(sD)0('=^ ) bit operati ons. The 
number of polynomials in <&' is s; their degrees are dominated by D and 
the bit length of their coefficients is dominated by aD'^^'^ ^ . 

2. If a rational point with coordinates of bit length dominated by i is returned 
at Steps ffiJ^T^ , the rational number computed at Step (7(j)ii ) has bit 
length dominated by £ + aD'^^'' ^ . 

Proof. From Theorem[2?71 Steps (|7all7bp are performed within as'~'^''^ '> D'^'-''^'> 
bit operations. The obtained quantifier-free formula is a disjunction of (sD)'^^'' ^ 
conjunctions. Thus the loop (Step l(7c|) ) makes at most (sZ?)^^*' ^ calls to SemiAI- 
gebraicSolve. Each conjunctions involves {sD)'^^'''> polynomials of degree 
in . . . ,Ak,B] with integers of bit length dominated by aD'^'^^ \ 

Thus, from Proposition [2TT1 Step (7(c)i| is performed within cr'^(^) (sD)'-'('^ ' 



bit operations and outputs a rational parametrization of degree D'^''^ ^ with 
integer coefficients of bit length dominated by aD'^^'^ \ Thus, the integers 
in the list computed at Step (fTdl) have bit length dominated by aD'^^'^ \ This 
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implies that the polynomial obtained from Steps l|7el|7gP has rational coefficients 
of bit length dominated by aD'^^^ \ Assertion ^ follows immediately. 

The bit complexity of these steps is obviously negligible compared to the 



cost of Step |7(c)i|. The substitution phase (Step FZhl) has a cost which is still 



dominated by the cost of Step |7(c)i|. As announced, the obtained formula $' 
contains s {k — l)-variate polynomials of degree D with integer coefficients of 
bit length dominated by aD'^'^^ ^ . □ 

We prove now Proposition 14.11 by induction on k. The initialization of the 
induction is immediate from Lemmata 13.51 and 14.31 

Suppose that k > 1. Suppose that the execution of FindRationalPoints(<i>) 
stops at Steps |[2l), or (0]) or ([6]). From Lemma ITSj we are done. Suppose now 
that we enter in Step ([7]). 

By Lemma rOlf l]). the formula $' computed at Step (|7hp contains s (fc — 1)- 
variate polynomials of degree D and coefficients of bit length dominated by 
(jDO(k ) g^j^j obtained within (t*^^^^ (sl?)'-'^'^ ^ bit operations. The induction 
assumption impHes that 

• Step |[7il requires cr'-''-^-' (sZ3)*^^'^^) bit operations, 

• If a rational point is contained in L (Step 1^), its coordinates have bit 
length dominated by aD'^'-'^ \ 



Hence, by Lemma H^H (2|) . the rational number computed at Step I7(j)ii 
bit length dominated by aD'^'^^ \ Moreover, the cost of Steps ( 7(j)ii||7(j)iii 
negligible compared to the cost of previous steps. 



has 
is 



5 Rational sums of squares 

Consider a polynomial / G Z[a;i, . . . ,a;„] of degree 2d whose coefficients have 
bit length bounded by t. If we choose v as the vector of all monomials in 
2.[xi, . . . ,Xn] of degree less than or equal to d, then we consider the set of 
real symmetric matrices M = of dimension D = ("jj^"^) for which / = 
v^ • M • V. By Gaussian elimination, it follows that there exists an integer 
k<\D{D + l)- ("+2'') such that 

M = {Mo + FiMi + . . . + YkMu, Fi, . . . , Ffc G M} (6) 

for some rational symmetric matrices Mq, . . . ,Mk. The polynomial / can be 
written as a sum of squares of polynomials if and only if the matr i x M can be 
completed as a symmetric positive semidefinte matrix (see Laurent ( 200lh ). Let 
Y = (Yi,...,Yfe), we define 

5 = {y G M'' I MiY) h 0, M{Y) = M{Yf, / = v^ • M{Y) ■ v}. (7) 

It is clear that S C M''' is a convex set defined by setting all polynomials in 

$(ri, . . . , Ifc) = {(-1)(*+^)to„ i^O,...,D-l} (8) 

to be nonnegative, where the m^'s are the coefficients of the characteristic poly- 
nomial of M{Y). The cardinality s of $ is bounded by D and <& contains poly- 
nomials of degree bounded by D whose coefficients have bit length bounded by 
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tD (see lPowers and WormanrJ 1 1998I )). Hence the semi-algebraic set defined by 
CD is 

5 = {(ri,...,rfe) eiR^- | (-i)(^+^)m, > o, o<i<d-i}. (9) 



The result below is obtained by applying Theorem ll.ll to the semi- algebraic 
set defined above. 

Corollary 5.1 Let f G Z[a::i, . . . , a;„] of degree 2d with integers of bit length 
bounded by t. By running the algorithm FindRationalPoints /or the semi- algebraic 
set defined in f^i, one can decide whether f is a sum of squares in Q[xi, . . . , Xn] 
within t'^'-^-'I?^''^ ^ bit operations. Suppose f = ft ■> fi ^ Qi^ii ■ ■ ■ i^n], then 
the bit lengths of rational coefficients of the fi 's are bounded by tD'^'^^ \ 

Remark 5.2 Applying Proposition \1.3\ bv Khachiyan and Porkolab to the semi- 
algebraic set defined in one can decide whether f is a sum of squares 
in Z[xi, . . . , Xn] within t'^(^) D'^'^'' ^ operations. Suppose that f — Yfi,fi & 
. . . , then the bit lengths of integer coefficients of fi are bounded by 

^£)0(fc*)_ 

Porkolab and Khachiyan showed that the non-emptiness of the convex set 
defined in over the reals can be determined in 0{kD^) + i:)0(min{fc,D"}) 
arithmetic operations over £j)0(™i {fc.^ ^^-bit numbers, wher e i is t he maximal 
bit length of the matrices Mi (see Porkolab and KhachivanI (1997)). Suppose 
S ^ i.e., / G Q[a;i, . . . ,a:„] is a sum of m squares in K[xi, . . . ,Xn] where K 
is an algebraic extension of Q. If is a totally real nun iber fie l d, the n / is also 
a surti of squares in Q[a:i, . . . ,a;„], i.e, ^flQ" fsee lmilad ll2009[ l: iKaltofenI 
(l2009fl1. The follo wing lemma and proof can be deduced from arguments given 



KaltofenI ()2009l ). 



Lemma 5.3 Suppose Q ~ (G, Gq, Gi, . . . , Gk) is a rational parametrization for 
the semi- algebraic set S defined in (Tj) computed by SemiAlgebraicSolve. Suppose 
d is a real root of G such that 

Y{^) = i(Gi(i?), G2(^), . . . , Gum e 5, (10) 

Then for any real root di of G, we have 

Y{§,) = i(Gi(i?,), G2{A), Gfc(i?,)) e S. (11) 

Moreover, i^ the polynomial G has only real roots, then the point defined by 
degG ^(^i) c- rational point in S. 

Proof. Since l"(i9) £ S, the matrix M(Y{'&)) is positive semidefinite. We 
can perform the Gaussian ehmination over Q(i9) to obtain the decomposition 
M(y(i?)) = A{§)'^A{^). It is clear that for any real root ^i of G, M{Y{^,)) = 
A{'&i)'^ A{'&i) is also positive semi-definite, i.e., Y{di) € S. Moreover, if G has 
only real roots t?i, then qj^.^^q ^A'^i) ^ 'Q- follows that the point defined 

dc-rG SS*^ ^(^i) is ^ rational point in S. 

eg ^ 

The above discussion leads to the following result. 
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Theorem 5.4 Suppose f € T\x\, . . . ,Xn]. There exists a function RationalTo- 
talRealSolve which either determines that f can not be written as sum of squares 
over the reals or returns a sum of squares representation of f over Q[xi, . . . , x„] 
if and only if the polynomial G outputted from the function SemiAlgebraicSolve 
has only real solutions. The coordinates of the rational coefficients of polynomi- 
als fi in f = ff have bit length dominated by tD'^^^^ and the bit complexity 
of RationalTotalRealSolve is t°^'^) D'^^''\ 
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